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Abstract: We analyze the effect of an isospin chemical potential /i/ in the Sakai-Sugimoto 
model, which is the string dual of a confining gauge theory related to large Nc QCD, 
at temperatures below the chiral symmetry restoration temperature. For small chemical 
potentials we show that the results agree with expectations from the low-energy chiral 
Lagrangian, and the charged pion condenses. When the chemical potential reaches a critical 
value /i/ = /icrit — 1.7mp, the lowest vector meson (the "rho meson") becomes massless, and 
it condenses (in addition to the pion condensate) for ^u/ > /Xcrit- This spontaneously breaks 
the rotational symmetry, as well as a residual U{1) flavor symmetry. We numerically 
construct the resulting new ground state for m > /Xcrit- 
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1. Introduction and summary of results 

Unraveling the phase diagram of QCD, in particular at non-zero baryon density, remains 
a major challenge (see [1] for a recent review). Although lattice simulations have provided 
us with a great deal of information about the finite-temperature behavior of QCD, it is 
technically much less clear how to reliably extract the physics at finite baryon chemical 
potential, because of the notorious sign problem. One way to get around this problem is to 
consider non-zero isospin chemical potential instead [2, 3]. Although there are no physical 
systems in nature in which only a large isospin density is achieved (because weak decays 
do not conserve isospin), studying such systems may nevertheless provide us with useful 
knowledge about the full phase diagram of QCD.^ 

The effects of a small isospin chemical potential m can be analyzed using the chiral 
Lagrangian [3]. For chemical potentials below the mass of the pion, one finds that the 

^Neutron stars have a large isospin density, but in this case the baryon density is also large, in contrast 
to the idealized situation studied here, in which the chemical potential and charge densities are turned on 
only in a U{1) subgroup of the SU(2) isospin group. 
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pion masses (as well as the masses of other isospin triplets) split: one of them increasing, one 
decreasing and one remaining constant. The ground state, however, remains unmodified 
from that of QCD at jii = 0, until m reaches the pion mass rriT^. At this point, one 
of the charged pions condenses, with a magnitude that grows with /x/. This situation 
persists until /xj becomes of the order of the mass of the first massive meson, where chiral 
perturbation theory breaks down, and it is not clear how to analyze what happens for 
larger values of (except by lattice simulations, but these are hard to perform for large 
pion condensates). One may consider toy models in which higher-derivative interactions 
or interactions with massive vectors are added, such as the Skyrme model, but it is not 
clear how to do this in a controllable way. Only at extremely large values of jii does one 
regain analytic control. Here, due to asymptotic freedom, perturbation theory applies and 
predicts a (u^^d) condensate, much like a fermion superfluid [4]. Given that the quantum 
numbers of the condensate are the same as those of the pion condensate for small /i/, it 
has been conjectured [4] that no phase transition occurs as /// is increased from 171-^ to 00. 
However, this remains a conjecture, and it has also been suggested that the rho meson 
may condense when the isospin chemical potential is of the order of the rho meson mass, 
spontaneously breaking the rotational symmetry [5, 6, 7] (see [8] for a review). So far, 
lattice studies are able to reproduce the results of the chiral Lagrangian (see, for instance, 
[9]), and studies arc currently under way to explore more of the intermediate and high 
isospin chemical potential regimes. 

In the present paper we analyze the behavior of QCD-like theories at moderate values 
of the isospin chemical potential, using an alternative analytic route. Instead of analyzing 
QCD, we will analyze the Sakai-Sugimoto model [10, 11] (see also the review [12])^. This 
model has a dimensionless parameter X^/R (we will explain the meaning of A5 and R in 
the next section). It is an easily- analyzed weakly-curved string theory in the limit of large 
A5/-R, while it goes over to large Nc QCD in the opposite limit; our analysis will be limited 
to the former limit, where the model contains many additional degrees of freedom beyond 
those of QCD. For any value of the dimensionless parameter, this model is dual to a large 
Nc confining SU (Nc) gauge theory with Nf massless flavors, which exhibits spontaneous 
chiral symmetry breaking. We will analyze the case of an isospin chemical potential, as 
opposed to a baryon chemical potential^, for three reasons. First, we will work in the large 
Nc limit, and in this limit the mesons (which are the relevant degrees of freedom for the 
isospin chemical potential analysis) are much more similar to the mesons of QCD than the 
baryons (which are relevant for the baryon chemical potential analysis); in particular they 
are weakly interacting and have finite masses. Note also that at large baryon chemical 
potential the phase diagram for large Nc is quite different than for Nc = S [19], but 
the phase diagrams seem similar for the isospin chemical potential case. Second, in the 
Sakai-Sugimoto model, mesons are much easier to control than baryons. The light mesons 
correspond simply to modes of massless fields living on the D8-branes, while baryons are 

^Flavor chemical potentials in other holographic models, which do not exhibit dynamical chiral symmetry 
breaking, were recently studied in [13]. 

baryon chemical potential in this model was recently analyzed in [14, 15, 16, 17, 18]. 
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complicated solitonic objects with a size of order the string scale [20, 21, 22, 23] (or maybe 
even smaller), so it is not clear how to write down the full effective action describing the 
baryons and their condensation. Third, as mentioned above, the phase diagram for a finite 
isospin chemical potential can potentially be compared with lattice simulations, because 
the "sign problem" which plagues simulations with finite baryon chemical potential is 
absent [2, 3]. 

The Sakai-Sugimoto model describes a chiral gauge theory with an SU{Nf) l x SU {Nf)ji 
global symmetry which is spontaneously broken at zero temperature. At finite temper- 
ature it exhibits (at zero chemical potential) three possible phases [24, 25]. In the low- 
temperature phase, the model is confined and the chiral symmetry is broken to its diagonal 
subgroup. In the intermediate-temperature phase, which only appears for sufficiently large 
constituent quark mass ^, the gluons deconfine but the chiral symmetry remains broken. 
Finally, in the high-temperature phase, the chiral symmetry is restored. The behavior 
of low-spin mesons, which we will focus on in this paper, is described by small fluctua- 
tions of the flavor D8-branes, which can be treated (for Nj -C Nc) as probes embedded 
in the D4-brane background. The meson spectra at zero chemical potential and non-zero 
temperature have been analyzed in [31]. 

Since the model has a vanishing quark mass, the pions are massless, and upon turning 
on an isospin chemical potential the charged pions immediately want to condense. In the 
absence of a pion mass, the pion condensate is stabilized by the non-linear interactions of 
the chiral Lagrangian; for any non-zero value of the isospin chemical potential, the pion 
goes to its maximal possible value, given that it is a Nambu-Goldstone boson which lives 
on a compact target space. The pion condensate leads to a small technical problem in 
our analysis; since the pion condensate in the Sakai-Sugimoto model corresponds to the 
holonomy of the D8-brane gauge field in the holographic z direction, one has to deal with 
solutions with Az 7^ 0, while usually it is more convenient to analyze this model in the 
Az = gauge. To avoid this complication, one can perform an SU{Nf)L x SU {Nf)ji global 
symmetry transformation (which is a "large gauge transformation" in the Sakai-Sugimoto 
model) which eliminates the pion vacuum expectation value, and then continue working 
in the Az = gauge. This global symmetry transformation maps a vectorial isospin 
chemical potential with a pion condensate into an axial isospin chemical potential with no 
pion condensate. Thus, instead of studying the effects of the vector-like isospin chemical 
potential in the presence of the pion condensate, we can equivalently study the effects of 
an axial isospin chemical potential in the trivial Az = vacuum. Using this method, we 
study the effects of the chemical potential as /x/ is increased. In particular, we determine 
the behavior of the masses of the low-spin mesons as a function of /i/. Note that axial 
flavor chemical potentials may be interesting in their own right for studying the full phase 

*The Sakai-Sugimoto model has another dimensionless parameter, which we will denote by L/i? in the 

next section, which determines the constituent quark mass (as measured, for instance, from the spectrum [26] 
or decay rates [27] of high-spin mesons). This parameter has a large effect on the physics in the limit we 
work in, though it is expected to decouple in the QCD limit. The bare quark mass (and, consequently, also 
the pion mass) is zero for any value of this parameter. Some recent attempts to add non-zero bare quark 
masses to this model appear in [28, 29], see also [30]. 
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diagram of QCD; as mentioned above, an axial isospin chemical potential is equivalent to 
a vector-like isospin chemical potential by a global symmetry transformation, but one can 
also study an axial chemical potential, and we will describe how the Sakai-Sugimoto 
model reacts to such a chemical potential as well. 

Our results for the Sakai-Sugimoto model at low temperatures (below the deconfine- 
ment temperature) and finite isospin chemical potential (in an SU{2) subgroup of the 
diagonal SU{Nf) subgroup of the flavor group) are the following. For small chemical po- 
tentials we correctly reproduce the results of the chiral Lagrangian, and the charged pion 
condenses to its maximal value. The isospin chemical potential m explicitly breaks the 
isospin symmetry to U{1); the pion condensate then breaks this U{1) symmetry, but in 
the vacuum with the pion condensate there is a new U{1) symmetry which is an axial 
subgroup of U {Nf )L x U{Nf )ii which is unbroken. In this phase the isospin charge density 
is exactly linear in /xj. The interesting physics starts to appear when the chemical poten- 
tial becomes of the order of the p meson mass. In this paper we analyze this physics for 
the special value of the parameters L = ttR, where the D8-branes are anti-podal to the 
anti-D8-brancs. We find that there is a critical value of m = /Ucrit ~ 1.7 where the p 
meson becomes unstable and condenses. This condensate breaks the rotational invariance 
from S0{3) to SO{2), and also breaks the remaining U{1) symmetry completely. We nu- 
merically construct the new ground state, and find that for m slightly larger than /Xcritj the 
condensate behaves as (p) oc y^^u/ — /Xcrit , as expected for a second order phase transition. 
On top of the p meson condensate, the pion condensate persists as well. The analysis of 
larger values of p,i, as well as of other values of L/R (where there is also a non-trivial 
intermediate temperature phase ^), is postponed to future work. In particular, it would be 
interesting to check the stability of the new state we find, and to see if additional mesons 
also condense as jU/ is further increased. In the Sakai-Sugimoto model (as in all weakly 
curved holographic models) the mesons of spin larger than one are much heavier (by a 
factor of a/ A5 /i?) than the low-spin mesons; however, when fij reaches the scale of the 
mass of these mesons, they may also want to condense. Our analysis is performed in the 
approximation of pj <C X^/R^; for larger values of pj one has to include also DBI and 
Chern-Simons corrections to the D8-brane effective action, and it would be interesting to 
see if these wash out the rho meson condensate that we found or not. 

Let us end with a few remarks concerning the relevance of our analysis for real (large 
Nc) QCD. The large isospin chemical potential analysis in QCD reviewed above, and the 
subsequent conjecture about the absence of phase transitions [4], crucially depends on the 
underlying asymptotic freedom of the theory. The Sakai-Sugimoto model, however, does 
not exhibit asymptotic freedom (except in its QCD limit), so apriori it is not clear how 



®The analysis of isospin chemical potentials in the high-temperature phase of the Sakai-Sugimoto model is 

a straightforward generalization of the analysis of [32, 33]; in this phase the dynamics involves a condensation 
of deconfined quarks, and there is no difference between the isospin and baryon chemical potentials, as well 
as between vector-like and axial-like chemical potentials. Isospin chemical potentials in the intermediate 
and high temperature phases of the Sakai-Sugimoto model were recently analyzed in [34], but the results 
there do not agree with ours. 
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this model behaves even for large ^u/. We do not know if the phase transition described 
in the previous paragraph persists when we go to the QCD limit of small X^/R or not. 
However, it is plausible [5, 6, 7] that a similar phase transition could occur also in QCD in 
the regime of intermediate m (of order the p meson mass), where the analysis of [4] does 
not apply. It would be interesting to compute corrections to our analysis in inverse powers 
of As/i? (related to a' corrections), to see if these tend to strengthen or weaken the rho 
meson condensate, as we go in the direction of approaching QCD. 

We begin in section 2 with a review of the Sakai-Sugimoto model, focussing on the 
realization of chiral symmetry and on the corresponding moduli space of vacua. In section 
3 we analyze the isotropic solutions with chemical potentials, and compare them (for small 
chemical potentials) to the expectations from the chiral Lagrangian. In section 4 wc analyze 
the stability of the isotropic solutions, and find that there is a critical value of the chemical 
potential where they become unstable; above that value we construct a new branch of 
solutions in which the rotational symmetry is spontaneously broken. An appendix contains 
a discussion of some technical issues in the numerical analysis. 

2. The configurations with no chemical potential 

2.1 Reviev\^ of the phases of the Sakai-Sugimoto model 

The Sakai-Sugimoto model [10] is the only known example of a holographic dual for a 
3 + 1 dimensional gauge theory with a continuous SU{Nf)£^ x SU{Nf)ji chiral symmetry 
which is spontaneously broken. It is derived by considering the decoupling limit of Nc D4- 
branes, compactified on a circle of radius R {x4 = X4 + 27ri?) with anti-periodic boundary 
conditions for the fermions [35], intersecting with Nj D8-brancs at X4 = and Nf anti- 
D8-branes at .T4 = L. On the field theory side this gives a 4 + 1 dimensional SU (Nc) 
maximally supersymmetric gauge theory (with a specific UV completion which will not be 
relevant for us), compactified on a circle with anti-periodic boundary conditions for the 
adjoint fermions, and coupled to Nf left-handed fermions in the fundamental representation 
of SU{Nc) localized at X4 = 0, and to Nf right-handed fermions in the fundamental 
representation localized at X4 = L. In this section we will briefly review the holographic 
description of this model (see [10, 24] for more details). 

The holographic dual description of this model involves (in the large Nc limit with 
fixed Nf) Nf probe D8-branes in a closed type IIA string background with a metric, 
Ramond-Ramond 4-form and string coupling given by 



the 4 -|- 1 dimensional gauge coupling by g'^ = {27r)'^gsls- The submanifold spanned by X4 
and u has the topology of a cigar with u > ua, and requiring that this has a non-singular 
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geometry gives a relation between ua and R, 



3 \ 1/2 



The background is weakly curved and can be well-approximated by gravity whenever the 
't Hooft coupling A5 = g'^Nc of the 4 + 1 dimensional gauge theory is much larger than R. 
In the rest of this paper we will work in this regime. Note that in this regime the mass 
gap is of the same order as the scale 1/R oi the Kaluza-Klein modes in the gauge theory, 
so the 3 + 1 dimensional theory does not decouple from the higher Kaluza-Klein modes (as 
it does in the opposite limit of A5 ^ R, which is equivalent to large Nc QCD at energies 
much smaller than 1/R). In this regime the confining string tension Ts ~ As/i?^ is much 
larger than the scale of the mass gap, enabling an approximation in which the higher string 
excitations are ignored. 

The probe D8-branes fill the 3 + 1 dimensional space-time and the four-sphere, and 
trace a curve u{x4) in the remaining two coordinates. This curve is a solution to the 
equations of motion of the probe D8-branes, with the boundary condition that u — > 00 
at X4 = (with one orientation of the D8-branes) and at X4 = L (with the opposite 
orientation). At 2:4 = L/2 the solution reaches its minimal value of u, which we denote by 
u = uq. In the special case of L = ttR, no = u\ and there is a simple solution in which u is 
independent of X4, with a brane at X4 = joining smoothly (at u = u\) with an anti-brane 
at X4 = ttR. In other cases the solution cannot be written down explicitly, but it satisfies 
the first order differential equation 

- UQ^/fiuo). (2.3) 



7(u) + 



3 „/2 



(Rm 
\ u J f{u 



As we will review below, the fact that the D8-branes and anti-D8-branes join together 
smoothly is a reflection of the spontaneous breaking of the (classical) U{Nf)L x U{Nf)ii 
chiral symmetry to the diagonal U{Nf)y. 

At finite temperature T we can describe the theory by a solution of the equations of 
motion in Euclidean space. There are two possibilities for the closed string background. 
For T < l/27ri2 the dominant solution is simply the Euclidean continuation of ( p.lD , with 
the time direction periodically identified t = t + 1/T. In this solution the X4 circle shrinks 
to zero size at u = ua while the Euclidean time circle never shrinks (see the left-hand 
side of figure |l|). For T > 1/2ttR the solution is the same but with the roles of t and X4 
interchanged, such that the f{u) factor in the metric sits in front of dt^. In this phase u\ 
is related to the periodicity in the Euclidean time direction, 

1 4Wi^l,4^^/^ (2.4) 



T 3 V ^^A 

In the low-temperature phase T < 1/2ttR the profile of the D8-branes remains exactly 
the same as at zero temperature. At higher temperatures there are two possible profiles 
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t X4 



t X4 



t X4 





confined phase 



deconfined with chiral 
symmetry breaking 



deconfined with chiral 
symmetry restored 



temperature 

Figure 1: The topologies of the background and of the D8-branes in the three phases of 
the Sakai-Sugimoto model. For L > 0.97R, the model jumps directly from the first phase 
to the third phase at T = \/2'kR. 



for the D8-branes [24, 25]. For 1/2ttR < T < 0.154/L (the "intermediate temperature 
phase"), the dominant profile is similar to the one above, and the branes join together and 
spontaneously break the chiral symmetry (see the middle of figure |l|). The function u^x^) 
obeys a very similar equation to the one above. 



u 



u, 



1 + 



Rda 



u 



7R 



(2.5) 



On the other hand, for T > 0.154/L and T > 1/2ttR (the "high temperature phase"), the 
dominant configuration is that of separate branes at 2:4 = and anti-branes at X4 = L, 
which smoothly end at the horizon u = u\ (see the right-hand side of figure |T|). In this 
phase the chiral symmetry is restored. 



2.2 The chiral symmetry and the moduli space 

In holographic dualities (generalizing the AdS/CFT correspondence [36]) there is a one- 
to-one mapping between gauge symmetries in the bulk and global symmetries of the dual 
field theory. More precisely, any gauge field ^^(x,u) (where x are the coordinates of the 
dual field theory, u is the radial direction and we assume that we have performed a Kaluza- 
Klein reduction over any additional compact dimensions) that is present near the boundary 
u ^ 00 corresponds to a generator of the global symmetry of the theory. The theory is 
invariant under gauge transformations 

^ + igdt^g-^ (2.6) 

with a gauge parameter g{x,u) (an element of the gauge group G) which goes to the identity 
element g = I at the boundary. "Large gauge transformations" in which g{x, u 00) = go 
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correspond to global symmetry transformations with parameter . If the vacuum is 
invariant under such a "large gauge transformation" then the global symmetry is unbroken, 
and otherwise it is broken. 

This correspondence applies both to gauge fields coming from closed strings and to 
gauge fields coming from open strings. In the D8-brane configurations we discussed above, 
there are two independent U{Nf) gauge fields living near the boundary of space-time; one 
gauge field (Aj^) living at X4 = and another (Aj^) living at X4 = L. Thus, these theories 
have a U{Nf)L x U{Nf)ji global symmetry''. Note that in this case the gauge fields near 
the boundary are only functions of 3 + 1 out of the 4 + 1 dimensions of the gauge theory, 
reflecting the fact that the global symmetry acts on fields which are localized in 3 + 1 
dimensions. 

In a configuration where the D8-branes are separate from the anti-D8-branes, as we 
have in the "high temperature phase", we can perform separate constant gauge transfor- 
mations (of the form ( |2.6| ) with g = go ^ U{Nf)) for the fields Aj^ and Aj^, and the theory 
is invariant under these transformations. Thus, in this phase the full global symmetry is 
unbroken. 

On the other hand, in a phase where the D8-branes and anti-D8-branes are connected, 
the fields Aj^ and Aj^ are limits of a single gauge field living on the D8-branes, and we 
cannot perform independent gauge transformations on each of them separately. We can still 
perform a constant gauge transformation g = go on the full gauge field A^{x, z) (here z is a 
single-valued coordinate on the D8-branes, such that z — > —00 as the D8-branes approach 
the boundary at X4 = 0, and 2; — > 00 as the D8-branes approaches the boundary at X4 = L; 
see § |3.3| for details) . Near the boundary this acts on both A^ and A^ by the same gauge 
transformation g^, so it corresponds to the vector-like global symmetry in the theory. If 
we start from a configuration of the D8-branes which has a vanishing gauge field A^ = 0, 
then this transformation leaves this configuration invariant, meaning that the vector-like 
symmetry is unbroken. 

However, suppose that we want to act with a non- vector-like global symmetry transfor- 
mation, in which the transformation on A^ goes to gi G U{Nf) near the boundary, while 
that of Aj} goes to gn G U{Nf) (with gji 7^ gi). We cannot achieve such a transformation 
with a gauge parameter g{x, z) which is a constant. The simplest possibility is to choose 
g = g{z), such that g is a smooth function with g{—oo) = gi and 5(00) = gR. However, 
such a gauge transformation does not leave our configuration invariant, since it generates 
Az = igdzg~^ which is non-vanishing. In particular, the gauge-invariant holonomy 



^If the u CO limit of g{x,u) gives a result which depends on x, then the transformation is not a 
symmetry at all. 

^The axial (7(1) symmetry is actually anomalous so it is not a global symmetry. However, in this paper 
we will only discuss the classical theory in the bulk (the leading order in 1/Nc) so this will not be relevant 
for us [37, 38, 10, 39]. 
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which was the identity matrix U = I before this transformation, becomes U = gj^ gn after 
this transformation. Thus, the theory is not invariant under global symmetry transforma- 
tions with gL 7^ (7_R, reflecting the spontaneous breaking of the chiral symmetry. 

Of course, even when the global symmetry is spontaneously broken, a global sym- 
metry transformation still takes us from one vacuum to another, equivalent, vacuum 
of the theory. The theory we are discussing has a (classical) moduli space of vacua 
{U{Nf) X U{Nf))/U{Nf), and we can go from any vacuum to any other vacuum by a 
global symmetry transformation. Usually one characterizes the vacua by a matrix U which 
transforms under the global symmetry as C/ — > gJ^^Ugn; we see that in the Sakai-Sugimoto 
model this matrix is precisely the holonomy ( p.T] ) we discussed above^. The fluctuations in 
U are usually written asU = exp {iTra{x)T"' / fj^) , where are the generators of U{Nf) and 
the "pions" Tra{x) are the Nambu-Goldstone bosons; in the Sakai-Sugimoto model these 
particles come from fluctuations of the lowest radial mode of Az (this statement depends 
on the gauge; for instance, in the U = I vacuum one could choose the gauge = and 
then the pions would come from the boundary conditions on the other components of 
instead [10, 11]). 

Since all the vacua of the theory are equivalent, it is usually enough to study the 
vacuum U = I which is described by the solution = 0. However, in this paper we will 
be interested in turning on deformations which explicitly break the global symmetry. Such 
deformations destroy the equivalence between the different vacua, so one should study their 
effect separately in each vacuum. Nevertheless, since all the vacua are related by the global 
symmetry, we can always study the effect of the deformation on a vacuum with some non- 
trivial U = gJ^^gR by performing the transformation back to the vacuum U = I, but taking 
into account the fact that the global symmetry transformation acts also on the deformation. 
Thus, instead of studying the effect of the deformation as a function of the vacuum, we 
can stay in the vacuum U = I and study the effect of the deformation as a function of a 
global symmetry transformation {gL,9R) acting on the deformation parameters. 



3. Isotropic solutions with flavor chemical potentials 

In this section we will add isospin chemical potentials to the model described in the previous 



section. We begin in §3.1 with a general discussion of chemical potentials in the Sakai- 
Sugimoto model. We continue in § |3.2| with an analysis of the effect of chemical potentials 
on the chiral Lagrangian, which is the low-energy limit of our model. In § p.3| and § p.4| we 
find isotropic solutions of the full Sakai-Sugimoto model with chemical potentials, and we 
show that they reproduce our expectations from the chiral Lagrangian. 

3.1 Flavor chemical potentials in the Sakai-Sugimoto model 

As alluded to in the previous section, a gauge field Af^ in the bulk of a holographic model is 
dual to a global symmetry current in the field theory. As for any other field in the bulk, 



*Note that in the other vacua the preserved symmetry is not the vector-hke subgroup of U{Nf) x U{Nf) 
but a different subgroup, defined by g^^Ugn = U. 
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the behavior of the gauge field near the boundary includes independent non-normalizable 
and normalizable terms; in the Sakai-Sugimoto model this behavior is given by (as u — > oo, 
in the gauge = 0) 



From here on we use conventions with = 0,1,2,3, i,j = 1,2,3. Again, as for any 
other field, the coefficient By[x) of the non-normalizable mode is interpreted as a source 
term in the Lagrangian of the field theory of the form Jd^x Bi,{x)J'^ {x), while Pv{x) is 
proportional to the vacuum expectation value of Ju{x). The freedom in performing gauge 
transformations in the bulk does not affect the source term in the Lagrangian due to 
the conservation equation d^J'^{x) = 0. Adding a chemical potential to the field theory 
corresponds to adding a constant source for J^, so it corresponds to boundary conditions 
of the form ( |3.1| ) with B^{x) = fi6^fl. 

As described above, the gauge field living on the D8-branes in the Sakai-Sugimoto 
model comes near the boundary at two independent positions, which we denoted z — > — oo 
and z oo. Thus, for this gauge field we have two independent chemical potentials 
Ajj{x,z — oo) = hlSu,o and A^{x,z — > oo) = URSufi, where and piR are matrices in 
the adjoint of U{Nf)i and U{Nf)ii, respectively^. 

In the high-temperature phase, we have two separate stacks of D8-branes. The stack 
at X4 = only feels fiL and reacts to it, and the stack at X4 = L only feels fiji and 
reacts to it. These chemical potentials lead to a non-trivial gauge field on the D8-branes 
which leads to a specific charge density poix) which is linear in the chemical potential 
(for small p), as described in [32, 33]. For small enough p (compared to AsT^) so that 
we can ignore the higher order terms in the DBI action, the solution in this phase has 
Fou,L = —3pLU^/'^/2u^/'^ and Fou,r = — 3//_R«^^^/2ii^/^, leading to a change in the action 



In this paper we will analyze the lower temperature phases, where the global symmetry 
is broken. It is then often convenient to change the notation from left and right chemical 
potentials to vector and axial chemical potentials, defined (in an arbitrary normalization) 
by PL = Pv — PA and pR = pv + pa- The vector- like chemical potential pv is even under 
z —z, while the axial-like chemical potential pA is odd. Note that the U{1) components 
of Pv and PA are the chemical potentials for the baryon number and the axial U{1), 
respectively. The non-Abelian components correspond to isospin-like symmetries, either 
vector or axial. When studying the isospin we will usually limit ourselves to an SU (2) 
subgroup of U{Nf), and turn on chemical potentials proportional to 173 in this subgroup; 
the generalization to arbitrary chemical potentials (in the Cartan subalgebra of SU{Nf)) 
is straightforward. 

Our definition of pv and pA is natural in the U = I vacuum where the vector-like 
subgroup is unbroken (see footnote ^). Under the U{Nf) x U{Nf) global symmetry, the 

^We will use conventions in which Aq is the Lorentzian gauge field, such that Ao and fi are real. 




(3.1) 



AS* oc — tr{p'j^ + pji)u 



(3.2) 
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chemical potentials transform as dZ^I^LgL and fiR — > g^HRQR- The U{1) compo- 

nents of the chemical potentials are invariant, but the SU{Nf) components are not. As 
discussed at the end of §2^, we can rotate any vacuum to f7 = / by an appropriate global 
symmetry transformation, but we have to remember to act with the same transformation 
also on the chemical potentials. In particular, studying the same chemical potential in 
different vacua translates into studying different chemical potentials in the U = I vacuum. 



3.2 Chemical potentials in the chiral Lagrangian 

The spectrum of mesonic states in the Sakai-Sugimoto model at zero temperature [10] 
includes massless pious (which are the Nambu-Goldstone bosons of the chiral symmetry 
breaking) and other massive mesons. At low enough energies the theory can be described 
just by the effective action of the pions; on general grounds this is given by the chiral 
Lagrangian, which depends on the matrix U which we described in the previous section 
(this was explicitly derived from the Sakai-Sugimoto model in [10])^''. For large enough 
temperatures one expects additional degrees of freedom to become important. However, in 
the low temperature and intermediate temperature phases it is still true [24] that the only 
light states which are charged under the flavor symmetry^^ are the pions, and there is a gap 
to any other charged states (in the low temperature phase the other states are just massive 
mesons, while in the intermediate temperature phase they include also massive deconfined 
quarks). Thus, in these phases that we are interested in, the low energy dynamics is 
still captured by the chiral Lagrangian. In particular, the response of the theory to small 
chemical potentials should also be captured by this Lagrangian. Thus, in this subsection 
we analyze the response of the chiral Lagrangian to chemical potentials. In the next 
subsection we will verify that the Sakai-Sugimoto model leads to the same results for small 
chemical potentials. Our discussion for vector-like chemical potentials is a special case of 
the discussion in [40] (with vanishing quark masses), but we will add also axial- like chemical 
potentials. As mentioned above (footnote since we are in the large Nc limit we will 
ignore the anomaly in the axial U{1). 

The chiral Lagrangian may be written as 

/:chirai = YTr(^.f/^"f/^)- (3.3) 

The field U is an Nf x Nj unitary matrix transforming under the global symmetry as 
U g'2^UgR. The covariant derivative D^, is just the normal derivative in the absence of 



At higher energies the effective theory of the pions includes also higher derivative terms. Including 
the leading higher correction, which comes just from the Yang-Mills action on the D8-branes, leads to the 
Skyrme model with a four-derivative term suppressed by the scale 1/ L. Note that this scale is of the same 
order as the mass of the massive mesons, so at this scale one cannot necessarily trust an effective action 
which includes only the pions. The full effective action of the Sakai-Sugimoto model includes also many 
additional higher derivative terms which are governed by the confining string scale, so at energies of the 
order of the string scale we can no longer use any effective action. 

^^As discussed in the introduction, we will not include the baryon number symmetry, which behaves very 
differently from the other flavor symmetries, in our analysis. 
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chemical potentials. However, in the presence of chemical potentials it is modified to 



Di,U = dyU - iSyfiiuLU -UfiR) = duU - i(5i.,o([/^y, U] - {^ua, U]), 
D,U^ = d,U^ + i6,,o{U^fiL - PrU^) = d,U^ - i5.,o([)Uv, U^] + {PA, U^])- (3.4) 

Note that the baryon number chemical potential fiy oc / does not appear in the action 
at all; this is not surprising since none of the light states in the model carry a baryon 
number charge. Thus, the baryon number chemical potential has no effect on the classical 
low-energy action. 

The potential energy arising from the Lagrangian (|3.3|) is 



Khirai = 4 Tr {{[pv, U] - {^A, f/})([/iy, U^] + {/UA, U^])] ■ (3.5) 



4 

The charge density is : 

PV = = - ^4^u{^lv + pa)u^ + u\^iv - pa)u -2^,v), 

PA = = Y^UifiA + Pv)U^ + U\fiA - t^vW + 2^a). (3.6) 

We can now analyze the dynamics in various cases. We can start by having a chemical 
potential just for the axial C/(l), ixa = fJ-al/^/^f- In this case the potential is simply a 
constant Khirai = ~fnPa SO there is no potential generated on the moduli space, and the 
axial charge density is given by pa = tr{pA)/ \/Nf = 2f^fia- The axial U{1) chemical po- 
tential has precisely the same, rather trivial, effect also when additional chemical potentials 
are turned on, so we will assume it is zero from here on. 

Next, let us turn on a vector-like isospin chemical potential /iy = p.ia^/2 [40]. Our 
potential then takes the form 

K^hirai = Tr{a^Ua^U^ - 1). (3.7) 

o 

This is maximized by configurations where U is along the identity and directions, namely 
f^max = e*"(cos(/3)I -|- i sin(/3)(T3 ) , and minimized by configurations in the ai^2 directions, 
f^min = e*"(cos(/3)(Ti -|- sin(/?)(T2). Thus, the minimal energy configurations in this case 
will involve matrices of the form C/min- This is natural since there are massless pions 
charged under the isospin symmetry, and they want to condense to their maximal possible 
value when we turn on the isospin chemical potential. The isospin charge density in the 
f^min configurations is given by tr(/0ycj3) = 2f^pi, while in the C/max configurations it 
vanishes. Note that in the vacua U = Umm the U{1) symmetry corresponding to the 
chemical potential that we turn on is broken. 

It is just as simple to analyze the case where we have both axial and vector isospin 
chemical potentials; assuming that both of them are in the us direction (/iy = pi(T^/2 and 
Pa = Pa,i(^3/2), we find a potential 

Khirai = ^[(//? - pXi) Tr {cT^Ua^U^) - (//? + plj) Tr (l)] . (3.8) 
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Thus, if fij > ^\ J the minimum is stih at the same C/min, and the vacuum energy turns 
out to be independent of fiA,i, while if jij < ijlj^ j then the minimum is at one of the 
f^max configurations, and the vacuum energy is independent of In all of these cases, 
the vector-like charge density vanishes in the U = I vacuum, while the axial-like charge 
density is tripAcrs) = 2f^fiAj- 

An important point is that, as we mentioned at the end of |2.2| , we can study all these 
vacua also by transforming them to the trivial vacuum U = I. Suppose we wish to study 
the theory with a purely vector-like isospin chemical potential, = I^R = fJ-if^s/'^- As 
discussed above, in this case U = I is a maximum of the scalar potential, so it corresponds 
to an unstable solution to the equations of motion. Xhe stable solutions are at C/ — fAnin- 
We can transform this to [/ = / in various ways. For example, we can choose gL = C^mim 
gji = I; this transforms the vacuum to U = I, and the chemical potentials to a purely axial 
form fj^L = —IJ-R = —ni(7s/2. Other choices give similar results (related by the U{Nf)y 
symmetry which is unbroken in the U = I vacuum). Thus, we can study the behavior at the 
minima of the scalar potential with a vector-like isospin chemical potential by analyzing the 
vacuum U = I with a purely axial isospin chemical potential, see figure ^. Note that this is 
consistent with the fact that [/ = / is a minimum of the potential in the latter case.^^ Note 
also that even though the pion condensate breaks the original U{1) symmetry preserved by 
a purely vector-like isospin chemical potential, there is a new U{1) global symmetry which 
is unbroken in each of the vacua U = Umin] after the global symmetry transformation to 
U = I this is simply the vector- like U{1) subgroup of the isospin symmetry (but in the 
original variables it is an axial- like symmetry). 

3.3 Isotropic solutions for small chemical potentials 

In this section we will look for solutions to the equations of motion of the Sakai-Sugimoto 
model in the presence of the various chemical potentials discussed above. As discussed in 



§3.1, such solutions involve configurations in which Aq is non- vanishing. In this section 
we will only look for solutions which are translationally and rotationally invariant. Such 
solutions depend only on the time t and on the radial coordinate z, and must have Ai = 0. 
We will start with the solutions for small fi, which can be compared to those of the 



chiral Lagrangian, and in § |3.4| we will generalize this to larger values of For small fi we 
can approximate the action for the gauge fields on the D8-brane by the Yang-Mills action; 
from the form of the solutions that we will find, it is easy to see that the DBI corrections 
to the YM action start becoming important when fj, ~ As/L^, so our solutions in this 
subsection are valid for /u <C X^/L? . Moreover, we will only work to the leading non-trivial 
order in /i, so we will ignore the backreaction of the gauge fields on the brane embedding 
u{x/i) which is of order /x^ (this is not present when L = ttR but is present in all other 
cases). Since we are interested in solutions with Ai = 0, the CS term will play no role, so 
we will ignore it in this section. 



^■^Note also that if we transform a general value of U (with a vector- like isospin chemical potential), which 
is neither of the form f/min nor of the form f/max, to U — I, we would get vector-like and axial-like chemical 
potentials that do not commute with each other; for such chemical potentials U — I is not an extremum of 
the scalar potential, so it is not a static solution to the equations of motion. 
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f/ = /, (tt) = 

IJ'R = Ulcc (73 (vector) 



f/ = /, (tt) = 

I^R = -fj^L oc ff3 (axial) 




Figure 2: Introducing a vectorial isospin chemical potential leads to an instability for 
two of the pions (top left), and the theory wants to sit at a new minimum U = Ua\in 
(bottom). A global symmetry transformation brings us back to the vacuum U = I; this 
rotation changes the vector chemical potential to an axial one (top right). 



As mentioned above, it is useful to choose a coordinate z which is single-valued along 
the D8-branes; we will choose to define such a coordinate by 



(3.9) 



where z (and du/dz) are positive for 2:4 > L/2 and negative for 2:4 < L/2 (z = corre- 
sponds to the tip of the branes at u = uq, while z = ztcxo correspond to the two asymptotic 
regions on the branes). Using this coordinate, we can write down the YM action density 
for the gauge fields on the D8-branes in the low temperature phase as 



5 low 
D8 



with T = ^TsR^j^'^uo{27ra')'^ / Qs (where, as in [24], Tg includes the tension of the D8-brane 
and the volume of the four-sphere) , and 

7(u) 



T 



dz 



R 



4^5/2 



8 km \/7. 



(3.10) 



u 



In the intermediate temperature phase, where the gluons are deconfined but the chiral 
symmetry is still broken, the action density reads 



Clint 
'-'D8 



T 



00 
dz 

—00 



4^5/2 



tr(F, 



+ ytr(Fo2: 



9^9/2 p 

'8N^ 17 



iT(Fl)+t,[Fl) 



1 

(3.12 
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Note that in the hmit uq = u\ these expressions simplify because then 7 = /^^ and 
/ = uqz^ /v? . These actions assume that the brane configuration is the same as it was in 
the absence of the chemical potential; as mentioned above, we expect the chemical potential 
to generate gauge fields of order which should back-react on the configuration of the 
brane at order /i^. 

As discussed above, we are looking for solutions with Ai = 0, so the relevant part of 
the action is given by the Fq^ terms, which appear with a measure 52 (-z) given by 



9 1- 
92{z) = 0-2^ X { 

i Uq 



\ul + UQzyl^jJj= lowT, 

' '^^ (3.13) 
(uq + Uq z^)^/^ \z\l/f'y intermediate T . 



The large \z\ behavior is obviously the same for both cases. 

For simplicity, we will work in the A^ = Q gauge. As discussed above, this means that 
we are at a specific point U = I m. the moduli space, but this point is related to all other 
points by the global symmetry transformation so we can still discuss the most general case. 
Naively, going to Az = Q gauge means that we are setting the full pion field to zero, and 
not just its zero modes, but in fact this is not true; as discussed in [10], in this gauge the 
derivatives of the pions show up through the boundary values of A^. 

Thus, we wish to solve the equations of motion with an ansatz of the form 

= , Ai = Q, 

(3.14) 

Aq = Al{z,t) + A^{z,t), Fo, = -d,Ao, 

where we separated the field Aq into its even part (under z —z) A^ (related to the 
vector- like symmetry) and its odd part Aq (related to the axial symmetry). As discussed 



in §3.1, for a solution with a given chemical potential, ^0 should approach as z — > —00 
and II R as z ^ 00, and the coefficient of the normalizable mode of ^0 a-t infinity (which we 
called pa{x) in ( |3.lD ), which is proportional to 1/z, is proportional to the charge density. 



With the ansatz ( |3.14| ), the equations of motion become 

aoa.Ao -i [^0,9,^0] = 0, 9,(52(^)^,^0) = 0. (3.15) 

The second equation implies that the general solution for the field strength takes the form 



92{z) Jo 52(5) \Jo 92{z') 

for some U{Nf)-valued functions fv{t) and fA(t), which must obey the additional equation 
of motion 

dofAit)-i[fvit)jAit)]=0. (3.17) 

Note that the contribution of fv to Aq is even under z — > — z (in fact it is independent 
of z) so it includes any vector-like chemical potential and charge density, while the con- 
tribution of /a is odd and includes the axial chemical potential and charge densities. In 
the special case of the low temperature phase with L = ttR such that uq = u\, it is easy 
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to perform the integral in ( |3.16| ) explicitly, and to write the general solution in the form 
Aq = 2/^(i) arctan(z/tiA)/vr + fv{t) (where one still has to impose ( |3.17] )). 

Near the boundary, the solution ( [3.16 ) approaches /y(t) ± fA{t) at z = ±00. In the 



case that the solution is time-independent, we thus identify fv = /Uy and /a = fJ-A, and 
we find such a solution whenever [/iv',/xyi] = 0. This agrees with the expectation described 



in §3.2, that for commuting fiy and jja, U = I should be a static (not necessarily stable) 
solution. By examining the terms going as 1/z near the boundary, we see that these 
solutions have pv = and pA = CpA for some constant C that we will compute below. 



Again, this precisely agrees with the expectations described in §3.2. Note that in the purely 
vector-like case the field strength vanishes identically. 

The conclusions of the previous paragraph hold both for the U{\) and for the SU{Nf) 
components; as in the chiral model, the U (1) components of the gauge field are independent 
of the SU{Nf) components so we can discuss them separately. In the SU{Nf) case, if we 
take Pv = Pif^s/'^ and pA = pA,icy^/2, then, as discussed in §^, the solution we find is 
either a maximum or a minimum, depending on the sign of //^ — p\ j. However, this does 
not affect the form of the solution, but just the mass spectrum of small fluctuations around 
the solution. 

As we also discussed in § |3.2| , when [pv,pa\ 7^ we do not expect to find a static 
solution, and indeed we see from ( |3.17| ) that there is no such solution, but that the time 
derivative of /a must be non-zero (even as z ^ ±00). This implies that the time-derivative 
of the pion field must be non-zero in this case, consistent with the fact that we are not 



sitting at an extremum of the pion potential. It is easy to write simple solutions to ( 3.17| ) 



which describe oscillations around the minima of the pion potential, but we expect that 
when we include interactions (at order 1/Nc) such solutions will eventually decay to the 
minima, so they are not very interesting and we will not discuss them further. 

In order to compute the precise coefficient C appearing in our expression above for the 
charge density, we need to compute the derivative of the action of our configurations with 



respect to pA (it is obvious that the action is independent of pv)- By plugging ( 3.16 ) into 
the action density, we obtain 

S = 2Tv{p\)( r ^) \ (3.18) 
so that the axial charge density is given by 



00 J ~ \ — 1 



Comparing this with our expressions of the previous section, we obtain 



92[z) 

This is precisely the same expression that one obtains by computing the kinetic term for 
the holonomy matrix U (defined in (|2.7D ) in the Sakai-Sugimoto model [10]; note that a 
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configuration with varying holonomy U is described in this model by 



F^.z = '^U-%U ( r . (3.21) 

52(2) \J-oo92[z)J 

For L = ttR we have [10] oc X^Nc/R"^- For L <^ R, on the other hand, the scale R 
drops out of the analysis and we have oc X^Nc/L^. Since it is not clear how to directly 
compute the chiral condensate in the Sakai-Sugimoto model, f-j^ is the best estimate that 
we have for the scale of chiral symmetry breaking in this model. 

We can now use the results above to analyze the response of the theory to various 
chemical potentials. In the case of a baryon number chemical potential, we see that its 
only effect is a shift in the constant mode of Aq, which has no effect on the physics. 
This is as we expect (for small fi), since we did not include in our analysis any fields 
which are charged under the baryon number symmetry. For vector-like isospin chemical 
potentials we saw that we have similar solutions with constant values of ^oj these solutions 
have no field strength, but they affect the theory through the commutator terms in the 
action. As discussed above, these solutions (with U = I) should be unstable maxima of the 
pion potential, and the spectrum of fluctuations around these solutions contains tachyonic 
modes. If we want to analyze the minima f/min of the pion potential with a vector-like 
isospin chemical potential then, as discussed above, we can transform them into the U = I 
configuration with an axial-like chemical potential. These configurations describe a charged 
pion which condensed to its maximal value (it is stabilized by the non-linear interactions 
in the chiral Lagrangian, associated with the finite moduli space of the pions). In these 
configurations there is a non-trivial field strength, leading both to a non-trivial free energy 
density and to a non-trivial charge density, coming from the condensation of the charged 
pions. These configurations are expected to be stable, at least for small chemical potentials 
where the truncation to the chiral Lagrangian is a good approximation; we will analyze 
their stability at chemical potentials of the order of the mass gap in section 

3.4 Isotropic solutions for large flavor chemical potentials 



The solution ( 3.16|) for axial chemical potentials found in the previous subsection is valid 
only for small fi, or more precisely, for fi <C As/L^. Our analysis ignored both the DBI 
corrections to the Yang-Mills action, and the fact that the back-reaction of the gauge fields 
which are generated will change the brane embedding. In this subsection we will generalize 
this analysis to larger values of /U. In general, the effective action on the D8-brane includes 
also additional corrections beyond the DBI action, involving derivatives of the gauge fields. 
In the configurations we will be discussing in this paper, the corrections coming from 
derivatives in the radial direction are small as long as the gravity approximation is valid, 
namely for A5 ^ R. Corrections coming from derivatives in space-time directions will be 
small as long as these derivatives are smaller than the string scale ~ X^/R^; at higher 
energies or momenta than this we must use the full string theory and we cannot use any 
effective action. However, for the static solutions described in this section, the DBI action 
is valid for arbitrarily large values of fj,. 
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We will solve the equations of motion here for a single D8-brane; as long as the chemical 
potential is a diagonal matrix, the full solutions (in the U = I vacuum) will be diagonal, 
so they will just be Nf copies of the single brane solution. 

In the low temperature phase, the relevant terms in the action for a single D8-brane 

are 

S ^ ^ ri.W^)Jl + i^)' - il^f"^ . (3.22) 

9s Jo \l \ u J f[uy 

For simplicity we will work here in the u coordinate instead of in the z coordinate that we 
used before; the translation between the two coordinate systems is straightforward. One 
equation of motion can be obtained from the "Hamiltonian" associated to translations 
in X4, 

\/ f(u) 

H = ^■'^ ' , (3.23) 

h I /^^D4y n^^-/(n)(Fo4)2 

V w y f{^? 

which is conserved, d^^H = 0. In addition, the equation of motion for the gauge potential 
Aq implies 

E = , ""^^"^ = const . (3.24) 



04 



Fo4 = |f-^) /(n). (3.25) 



Dividing (|3.24| ) by (|3.23|) we get for the electric field 

H \iiD4 

Note that i*04 is symmetric under z —z, which means that Fqz is anti-symmetric, so 
this is an axial-like solution Of course, the vector- like solution is still given simply by 
Aq = const, E = 0. 

In order to express this in terms of our field theory parameters we need to find the 
shape u{x4). The brane minimum is by definition at u = uq, where u' = 0. This fixes H 
in terms of uq to be 

H = ^ulf{uo) + E^f{uo) (-^y . (3.26) 
We can then solve for u' from the expression for H, which yields 



I / \3/2 . 



(3.27) 



^^This axial solution was also obtained in [32], but incorrectly interpreted as a vector-like baryon chemical 
potential. When the configuration is made symmetric under z —> —z, the resulting gauge field is no longer 



continuous at the tip of the brane [33] . 
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As before, for a vector-like chemical potential the solution is simply Aq = jiy and 
-^04 = 0. For an axial- like chemical potential^^ we have the following relations between the 



parameters of the field theory, and L, and the parameters E and uq of the solutions 
found above : 



1 



dxidiA^ 



L 



-E 



dxA 



du 



E' 



21 ^ 



Mo 



2H 



du 



(RdA 

fin) \ u J 



3/2 



/?3 



E' 



-1/2 



(3.28) 



-n.D4 



-1/2 



(3.29) 

where in the last equation H should be expressed in terms of uq using (|3.26 ). It is easy to 
see that by expanding these solutions to leading order in /i (and in E) one reobtains the 



solutions we discussed in the previous subsection ( 3.16 ). 

Note that if we have generic isospin chemical potentials with Nj > 2 and the values 
of iJ,\ are not equal along the diagonal, then generally the D8-branes will no longer be 
overlapping but will have different functions u{x4) except near the boundary. This was not 
visible in the previous subsection since the separation starts at order jj.^ for small /x. The 
only case where this separation does not occur is for L = ttR (i.e. for uq = u\), since then 
constant values of X4 always solve the equations of motion. 

The analysis in the intermediate temperature phase is similar. The action is now 



/2 _ z?2 ^ 
^04) ■ 



The Hamiltonian takes the form 

H = 



(3.30) 



(3.31) 



and the conserved E takes the form 



E 



u 



"04 




i?D4 



U 



12 



const . 



(3.32) 



Dividing (|]3|) by (|3T|) we get for the electric field the same expression (|3.25| ) as before, 
and the expression for H{uq) is also identical to the previous expression ( 3.26]) . Solving for 



^"^We remind the reader that this case is equivalent to a vector-like isospin chemical potential when the 
theory is expanded around the minimum of the pion potential. 



- 19 - 



u' we now find 



L^f(u) - ulf(uo) + {f(u)u^ - f{uo)ul) (^-^) ^ ■ (3.33) 



The relations between the parameters now take the form : 

3/2 



^D4 



-1/2 



(3.34) 

-1/2 



(3.35) 



4. Instability of isotropic solutions and rho meson condensation 
4.1 The fluctuation equations 

In the previous section we found the static isotropic solutions of the Sakai-Sugimoto model 
in the presence of flavor chemical potentials. These solutions are the unique isotropic 
solutions in the absence of sources, but there may exist additional non-isotropic solutions, 
which could dominate the thermodynamics for some values of the chemical potential. We 
will focus on the isospin chemical potential case. Since the model contains massive particles 
charged under the isospin symmetry, one might expect that for a large enough chemical 
potential these particles will want to condense and modify the solution, and we will see 
that this is indeed the case. 

In order to analyze the stability in general we need to use the full DBI action, and 
include also the Chern-Simons term in the action of the five dimensional gauge field (which 
includes quadratic terms in the fluctuations around the solutions of the previous section). 
However, as discussed above, for /x ^ ^b/L"^ the DBI action can be approximated just by 
the Yang-Mills action, and in the same regime one can also show that the contribution of 
the Chern-Simons term is much smaller than that of the Yang-Mills terms, so it can also 
be ignored. Thus, in this section we will use the Yang-Mills action to look for instabilities, 
and we will see that the isotropic phase discussed in the previous section becomes unstable 
already for /x which obeys jj, <^ \^ / L? ^ so that the analysis is self-consistent. For simplicity, 
we will consider in this section only the low-temperature phase. We will also consider only 
the special case of L = -kR for which the equations (and the isotropic solution) simplify. It 
should be straightforward to generalize our analysis, but we postpone this to future work. 

To analyze the stability we need to consider the spectrum of fluctuations of all the 
flelds in the theory around our solutions. In the limit of weak curvature, we can limit our 
analysis to the massless flelds living on the D8-branes, since the other mesons are much 
heavier. Because the chemical potential breaks Lorcntz symmetry, there are four types 
of fields which have to be discussed: there are KK towers of scalars, transverse vectors, 
and longitudinal vectors, and in addition there are the pions. All of these except the 
scalars arise from the five-dimensional vector field. In order to keep the vector equations 
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of motion reasonably compact, we will temporarily revert back to the u coordinate system. 
The action (in the approximation described above) then reads 



(4.1) 

As in the previous section, we will analyze the effect of a vectorial isospin density in 
the U = f/min vacuum where the pions have condensed by performing a global symmetry 
transformation taking this back to the U = I vacuum, in which the isospin chemical 
potential becomes an axial isospin chemical potential. So, we will analyze the stability of 
the solution constructed in § |3.3| (which we write here just for the D8-brane branch z > 0), 



(4.2) 

with J\f chosen to make limu_^oo ^{u) = 1- We will continue to use the gauge Au = 0. In 
this gauge the expansion of the gauge field in terms of fluctuations is 



Ao = A^{u)+5Af\oj,k,u)aae'^'+'^-\ A, = ^(u;, ^, n) e*-*+*^■^ A, = . (4.3) 



The transverse fields coming from 5Ai{L0,k,u) will be degenerate since there is still rota- 
tional invariance, but they need no longer be degenerate with 5Aq. 

The equations of motion for the diagonal gauge field fluctuations SA^^^''^'* and SA^q''^^ 
do not change upon inclusion of the chemical potential and the background ( |4.2D ; their 
spectra remain unmodified (in the approximation of /i/ <^ As/L^ where we ignore the 
Chern-Simons term in the action). For the other matrix elements, the action ( |4.lD leads 
to the following equations for the spatial components. 



+ 2^lIA{u)^o5Af^ - ^iA{u)di6A^i^^ = -^l]A{uf5Af'^ , (4.4) 
- 2iiiA{u)dQ5Af'^ + ^iIA{u)^i5A^^^ = -ii]A{uf8Af^ . (4.5) 



Here = d^d^ uP' — k"^ and F^y = d^Ai, — d^A^ — i[A^, Ay], so that with our conven- 
tion A^ = A^^'^aa we have F^^f! = d^A^'f' - d^A^^'^ + 2€''^^ A^ A'f^ . In addition, we have the 
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equations for the time- like and u-components of the gauge field, of the form 

+ {did'5A\^^ - do{d^5Af^) + ^iid^dAf^ A{u)^ = , (4.6) 
+ (did'SAf - do{d^5Af^) - ^iid^dAf'^ A{u)] = , (4.7) 



^0 

duidodA'i^ - di5Af^) - ^lIA{u)^u5A'i^^ + ^lI5A'i^A'{u) = , (4.8) 
a,(ao54') - di5Af^) + ^,iA{u)du5A'^^^ - ^,M^^^{u) = . (4.9) 

The equations of motion for the scalar fields will be discussed in the next subsection. 

Naively, one may expect that for positive /// the masses of the fields carrying posi- 
tive isospin charge would go down and the masses of the fields carrying negative isospin 
charge would go up (at least for small /i/). The equations above are indeed symmetric 
under m — A*/, A^'^'^ — > —A^'^'> which changes the sign of the isospin chemical potential 
and also exchanges the positively and negatively charged mesons. However, there is an 
extra symmetry of the equations (and of the solution (^)) under nj —fij together with 
z — > — z (when the equations are written in the z coordinate), which implies that in fact the 
spectrum is the same for positive and for negative (note that together with the previous 
symmetry, this implies that the spectrum is the same for the positively charged mesons as 
for the negatively charged mesons). The reason for this is the condensation of the pions; 
the expectation above holds when we expand around the trivial vacuum U = I, but it 
does not hold when we expand around the vacuum in which the pions have condensed (or, 
equivalently, when we expand around the vacuum U = I with an axial isospin chemical 
potential). Note in particular that the mesons that we analyze in this section carry a C/(l) 
charge in the new vacuum, but this is not the same as the original U{1) isospin charge to 
which the chemical potential was coupled^^. In any case, we will see that the masses of 
some of the fields (of both positive and negative U{1) charge) go down with fij, and that 
eventually they condense. Because of the symmetries discussed above, it is enough to con- 
sider positive values of fij, and to consider either the positively charged or the negatively 
charged mesons (but not both), and we will use this when solving the equations. 

4.2 Transverse vectors and scalars 

The analysis of the transverse vectors and the scalars is similar to their analysis without the 
chemical potential [10]. Let us first focus on the solutions of the 6 A- ' fluctuations, which 
describe massive transverse vector mesons. For these modes, it is consistent to set 6Aq = 
and then also impose the standard constraint for massive transverse vectors di5A^ = 0. This 
automatically takes care of equations (|4.6|)-(|49|). The spectrum of eigenfunctions A^^'^^ 



^^In the original variables, if we choose Umin oc ai, the equations above involve the ai and as components 
of the matrices. 
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Figure 3: The frequencies of zero spatial momentum solutions for the transverse vec- 
tors (left) and scalars (right), as a function of the isospin chemical potential, in units 
of \/ua/R'^4^ = 2/(3i?). The frequencies at /n/ = reproduce the spectrum found in [10]. 



of ( [4. 4]) and (4^) and the corresponding eigenvalues uj can be found numerically, using 



standard shooting methods, as a function of /x/. For the special case of L = ttR they 



are displayed on the left-hand side of figure in units of yu\/R^^ = 2/{3R). We see 
that the frequency of the lowest transverse vector goes to zero at a chemical potential of 
magnitude /icrit ~ 2.1 * 2/{3R) 1.7 rup, indicating a potential instability there. 

A similar story to the one above holds for the scalar mesons, arising from the scalar 
field y on the D8-brane^^. The scalar mesons feel the chemical potential through the 
commutators in the covariant derivatives. The equation of motion for the scalars in the / 
and (T3 directions are unaffected by the chemical potential (f4.2D , while the equations for 
the (Ti and (72 components are easily generalized from those in [10]. In the z-coordinate 
system they read (for L = ttR) 

where K{z) = (u/ua)^, and where the plus/minus sign refers to the y^^^ and y^^^ functions, 
respectively. These equations are diagonal in the basis of positively and negatively charged 
mesons, given by 

The resulting behavior of the uj eigenvalues as a function of chemical potential is displayed 
on the right-hand side of figure ^. We observe that the scalars also go to zero frequency 
leading to potential instabilities, but only for chemical potentials which are larger in magni- 
tude than those for which the first vector meson becomes unstable, so the vector instability 
must be analyzed first. 



^® There are additional scalar fields coming from the components of the gauge field along the S"*, but 
they are heavier in the absence of the chemical potential, so we will assume that they do not condense and 
ignore them in our analysis. 
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4.3 Pions and longitudinal vectors 



Let us now turn our attention to the pions and longitudinal vectors. For any non-zero 
momentum, these particles carry the same quantum numbers under the rotation symmetry 
preserved by the chemical potential, so it is possible for them to mix as the chemical 
potential is increased. Both the pions and the longitudinal vector excitations involve the 
two components 



Af'"^^ = ikiA'r^'^' , and 



,(1,2) 



(1,2) 



Again, the equations are diagonal when we choose 



6A^ 



(1) 



±i6A, 



(2) 



(4.12) 



(4.13) 



(2) 

We will use the polarization corresponding to the plus signs below, and keep 6Aq and 



SA^^'^ as the two independent components (with the others given by ( 4.13| )); as discussed 



above, the other polarization gives identical results. 

The boundary conditions determine how the pions and longitudinal vectors are encoded 
in these fields. The general solution to the equations of motion of the two functions At 
and ^0 near the boundary takes the form 



SA. 



(1) 



Cl 

C0 + - + 

z 



5A, 



(2) 



do + — + 

z 



(4.14) 



Here cq and do are the coefficients of the non-normalizable modes, while ci and di are 
the coefficients of the normalizable modes. The equation of motion ( [4.8| ) imposes one 
constraint on these four free parameters, thereby leaving a two-parameter family after the 
overall normalization is fixed. 

The solutions for the pions are expected to have a non-zero value of 5At at z — > ±00 
(as for vanishing chemical potential in the A^ = gauge), and we normalize this value to 
tt/2. On the other hand, the gauge field FQi must vanish as 2; — > ±00, for the solution to 
have finite action. This constraint determines in terms of cq, leaving 



pion 



6A!f\z 



00 



6A^^\z 00) 



TT Cl 

TT Z' Cl /l^ 

2 \uj + fij 



(4.15) 



+ 



Note that as — > we have SA^^^ = uo6A^\ as expected due to the enhanced Lorentz 
symmetry there. 

On the other hand, the solutions corresponding to longitudinal vector mesons are 
characterized by the absence of non-normalizable components in Ai and Aq, so that we 
have (in an arbitrary normalization) 



longitudinal vectors : 



5A^^\z^ 00) = i + , 



5A'^\z^ 00) 



1 



(4.16) 



UJ + flj z 



+ ... . 
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Figure 4: The energy w as a function of the isospin chemical potential /x/, in units of 
y^u\/R'^^ — 2/{'iR), for the sector consisting (at ^/ = 0) of the pion tt and the first two 
longitudinal vector mesons pL and ai. At /.t/ k, 2.1 * 2/(3i?) the first longitudinal vector 
meson becomes unstable. This coincides with the instability of the transverse vector 
meson depicted in figure |^. 



Note that if ci — > oo in ( [4.15| ), the pion modes look like longitudinal vector modes (with 
infinite overall normalization); this will be important for us later. 

Let us now discuss the relation between the boundary conditions at z = ±00. At /i/ = 
the pion is encoded in A^^ which is anti-symmetric in z; the relation of the pion to a chiral 
symmetry transformation implies that this should remain true near the boundary also for 
non-zero m, 

A^^\z = +00) = -A^^\z = -00) . (4.17) 

For the zeroth component, the boundary condition comes from requiring FQi = at the 
boundary, which gives 

^[,^^ {z = +00) = {uj + m) Ai^^ (z = +00) , 

(4.18) 

a''^\z = -00) = {uj - ni)Ai^\z = -00) = {-Lu + ij.j)A^^\z = +00) . 

Before we discuss the states at generic values of the momentum k, let us first analyze 
the k = case. This case can in fact be solved analytically, again specializing to the case 
of L = ttR, where A{z) = 2arctan(2;/iiA)/vr. The equation of motion (4^) becomes 

-2ni5A^Q^ + ua{1 + z'^ /uI)(^{ttuj + 2ni arctan{z/uA))dz6A''Q\z)^ =0. (4.19) 

This is solved by 

mI^^ = {uj + — arctan(z/nA)) • const. (4.20) 
vr 

In order to make this satisfy the boundary conditions (|4.18| ) we are forced to choose oj = ^j. 
This simple k = state becomes the pion in the limit iv, fii — > 0. We see that its frequency 
goes up with l/x/l and never vanishes. 
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In order to analyze the fate of the remainder of the spectrum in the pion/longitudinal 
vector sector, we resort to a numerical solution of the equations of motion. Modes satis- 
fying ( [4. 161 ) are found in the usual way with a one-parameter shooting method, starting 
from z = oo and looking for values of uj such that the solution goes to zero as z — > — oo. 
For the modes which are pions at /i/ = 0, i.e. modes satisfying ( [4.15 ), a search through the 



two-dimensional space spanned by ci and lo is required. Some details about this procedure 
can be found in the appendix. 

The upshot of this analysis is the following. The modes which start out as longitudinal 
vectors when fij = all have decreasing uj as /// is increased, and the frequency of the 
lowest one goes to zero at fij = ^crit- For the modes starting out with pion boundary 
conditions, the situation is more subtle. They start by going up, but they mix with the 
other modes. The low- momentum modes starting out as pions develop a large value for ci, 
which goes to infinity as /U/ /^crit- These modes therefore go over from pure pions to 
pure longitudinal vectors^^ as /i/ goes from zero to ficvit- Similarly, the longitudinal vectors 
develop a pion component, and as fii /Ucrit one linear combination of the modes which 
starts out as longitudinal vectors may be identified with the pion. 

In figure |^ we plot the value of uj for the pion and longitudinal vector states for a range 
of values of at fixed k. This plot first of all shows that a small chemical potential leads 
to massive pions with uj'^ = k'^ + fi'j, in accordance with our analysis of the chiral Lagrangian 
in the previous section. As is increased, the frequencies of the zero-momentum pions 
and of the lightest longitudinal vectors cross each other, and the lightest state becomes 
a longitudinal vector rather than a pion. Note that for non-zero k the two states cannot 
cross each other since they then have the same quantum numbers, but instead they mix as 
described above, and they come closer and closer to level-crossing as k ^ 0. 

For a value of m ~ 2.1 * 2/(3i?), the lowest-mass longitudinal vector meson (which we 
call the "rho") goes to zero frequency (at zero spatial momentum), indicating a potential 
instability. This is precisely the point where we also found that the lowest-mass transverse 
vector meson becomes unstable. Thus, the conclusion is that an instability of the rho meson 
develops at ficrit ~ 2.1 * 2/{3R) ~ 1.7 mp, where nip is the mass of the rho meson at zero 
chemical potential. The fact that the transverse vectors and longitudinal vectors go to zero 
frequency at the same value of fij is not surprising, since for zero spatial momentum (where 
we find the first instability) there is no distinction between transverse and longitudinal 
vectors. 

Note that the Skyrme model, which (as mentioned above) is a good approximation to 
the low-energy physics of the Sakai-Sugimoto model, indicates that the first instability (the 
first mode to go to zero frequency) as we increase /i/ comes from pions with /u/ and k of the 
order of the scale of the Skyrme term (which is of the same order as the rho meson mass 
in the Sakai-Sugimoto model). This differs from our conclusion that the first instability 
comes from vector mesons of zero momentum. Of course, this is not a contradiction, since 
the Skyrme approximation breaks down when fij and/or k become of order the mass of 



^'^Recall that to study the pions we arbitrarily fixed co = tv/2, but a more appropriate normalization 
when ci — > oo is Co ^ with ci remaining finite. 
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the vector mesons, and the additional vector fields must be included in the analysis then 
(as we did in our analysis). 



4.4 Vector meson condensation 

We have seen that at fij « 1.7 m^, the zero momentum modes of the "rho" vector meson go 
to zero frequency. This suggests that the theory becomes unstable there towards forming a 
condensate of rho mesons, in addition to the pion condensate which is already present; this 
was conjectured to happen also in QCD [5, 6, 7]. This condensate is of a vector field, so 
it breaks the rotational symmetry SO (3) SO (2), and since it involves the o"i and/or a2 
components of the meson, it also completely breaks the remaining U{1) flavor symmetry.^^ 
In order to find this new ground state, we assume without loss of generality that the 
third component of the vector rho meson condenses, and that the condensate is proportional 
to O"!. We also assume a static and translationally-invariant solution. In this convention, the 
Goldstone bosons for the broken rotational symmetry are the fluctuations of Ag^^ and A^^ , 

(2) 

while the fluctuations of ' include the Goldstone boson for the broken flavor U{1). We 
then have to solve the non-linear equations of motion 



3 J ^0 

(4.21) 



with the boundary conditions Af^ {z = ±00) = ±////2, A>^ {z = ±00) = 0. We expect 
that for ^[ < //crit these equations admit only the solution of the previous section which 
describes the pion condensate, i.e. 

/i/ < /icrit : ^n^^ = — arctan(— ) , Ai}^ = . (4.22) 

(3) 

For hj > /Xcrit we expect the appearance of an additional solution, which has both Aq 
and A^^ non-zero. 

We again resort to a numerical analysis in order to find the solution. Details of this 
procedure can be found in the appendix. The upshot is that indeed a rho meson condensate 
forms for fij > /Xcrit, so there is a second order phase transition there. Around fij = /icrit 
this condensate is, to good approximation, given by 

for fli < /Xcrit, 

(P) « <; (4.23) 

- /icrit for /X/ > /icrit- 

This condensate is drawn in figure |5|, using an arbitrary normalization in which we read it 
off from the coefficient of the normalizable mode of ^3^'*; it should be possible to translate 




^**There is another possible condensate of the vector meson, in which the unbroken U(l) is a diagonal 
subgroup of the U(l) flavor symmetry and an SO(2) rotational symmetry [7]. Without loss of generality, 
we can assume that in this phase the non-zero components of the gauge field are A^^' = A^^ (in addition to 
A^^). We have verified numerically that in our model this alternative condensate always has higher energy 
than the condensate discussed below, so we will not consider it further. 



-27- 



Figure 5: Appearance of the new ground state in which the p meson has condensed. 
This new ground state also still has a non-zero pion condensate. We plot (— C3) (in units 
of ma) and (—da) (in units of (ua/Rda)'^^^) for both solutions, where C3 and ds are defined 
by = ^lI{l/2 + cs/ttz + •••), 4'^ = ds/z + ■■■ ■ 

this to a physical normalization where the meson has a canonical kinetic term, but we 
have not done this here, and we do not expect this to make a big difference near the phase 
transition. The value of (— /i/ * C3) in this figure is the isospin charge density (up to a 
constant); note that in the solution with the rho condensate it rises more rapidly than in 
the solution that has only the pion condensate. 

The fact that the expectation value of the rho meson scales as ^/i/ — /^crit near the 
transition is expected on general grounds, by a Landau-Ginzburg type analysis of the rho 
meson effective theory, assuming that the effective mass squared of the rho mesons near 
the transition is proportional to (//crit — A*/), and that there is an effective potential of 
fourth order in the rho mesons which stabilizes the condensate. Our analysis implies that 
the coefficient of this fourth order term is positive. 

4.5 Potential generalizations 

As indicated above, our stability analysis so far is limited to the special case of L = ttR, 
and it would be interesting to generalize it to other values of L/R, and to see whether the 
instability of the isotropic phase persists for all values of L/R, and how //crit depends on 
this parameter. For small enough values oi L/R there is also an intermediate temperature 
phase in which the branes are connected [24], and it would be interesting to see if a similar 
instability arises in that phase as well or not. In that phase there are also deconfined quarks 
in the spectrum, which may play an interesting role. 

Our analysis so far is limited to the first instability which arises as we increase 
which leads to a condensation of the lowest vector meson. It would be interesting to see 
what happens as is increased further - do additional fields also want to condense? 
This analysis is more difficult than the analysis we performed here since the solutions for 
fJ-l > Z^crit which one has to expand around are only known numerically, but still it should 
be possible to do it. Of course, eventually, as we increase /Li/ to order As/L^, we also need 
to include the DBI and Chern-Simons terms in the action. Baryons may also eventually 
become relevant for the analysis (recall that in the Sakai-Sugimoto model, the ratio of 
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the baryon mass to its isospin scales as X^/R compared to that of the mesons [10, 21], so 
baryons are not expected to be relevant for small It would also be interesting to look 
for possible non-translationally- invariant solutions. 

are never 



Naively, our analysis in this section implies that the large /i/ solutions of §3.4 



relevant, since the isotropic solutions found in that section do not dominate already for 
smaller values of ^i. However, there are several situations where these solutions may still 
be relevant. One is the case of an axial U{1) chemical potential (which is also part of the 
full phase diagram of QCD, and may be interesting in particular at large Nc where the 
axial symmetry is approximately conserved). In this case there are no instabilities for small 
values of /i (for which the spectrum is unmodified), but the Chern-Simons term may well 
lead to instabilities at larger values of /i (of order \^/ L?), similar to those discussed in [41]. 
To see this one has to expand the action of a single D8-brane around the solutions of § [3.4| , 
using the full DBI plus Chern-Simons action. Another case where these solutions could be 
relevant is the case of Nj > 2 with generic SU (Nf) isospin chemical potentials, with all 
eigenvalues different. When L < ttR this leads to a separation of the D8-branes, which 
gives a mass to the off-diagonal matrix elements, and this may eliminate the instabilities 
discussed in this section which come from off-diagonal elements of the meson field. Again, 
new instabilities may then arise for fj, ~ As/L^; it would be interesting to investigate this 
further. 
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5. Appendix: physical modes 



In this appendix we provide some more technical details related to the determination of 
the meson spectrum, as well as the construction of the new ground state in which the 
rho meson is condensed. This analysis is, in both cases, essentially based on a standard 
shooting method for the two-point boundary value problem. However, there are a few 
interesting features which we would like to highlight. 

Let us first discuss the construction of solutions with asymptotic expansion (4.15) and 
boundary conditions ( 4.1S| ). When solving for the values of ci and uj for which such a 
mode exists, by shooting from z = oo to z = — oo, it is instructive to plot two curves, one 
for which the Aq^^ boundary condition is satisfied at z = —oo and one for which the A^rp 
boundary condition at z = — oo holds true. This leads to figure |6|. The intersection points 
of these two types of curves indicate physical states. The intersection which is almost on 
the w-axis has a 1/z coefficient (in A^) which is of similar magnitude as the asymptotic 
value of A^ (and thus corresponds to the pion). The other intersection point, away from 
the w-axis, is a longitudinal vector (when the 1/z coefficient of these modes is normalized 
to one, their value at z = ±00 goes to zero as A; — > 0). Following the figures, there is a 
clear sense in which a state which starts out as a pion eventually turns into a longitudinal 
vector, and vice versa, as fij is increased. As discussed in § |4.3| , this transition is smooth 
for any non-zero momentum. 
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Figure 6: The location of the physical modes for small and fixed k [k — O.Oby^UA/ R^^) 
for various values of /ij (in units of \fu\/B^ — 2/(3i?)). The horizontal and vertical 
axes label u (in the same units) and the coefficient of the 1/z term in (in units of u\, 



in the normalization of (4.15)), respectively. The blue, almost straight curves indicate 



where the boundary condition on A^^z = —00) is satisfied, while the red curves do the 



same for A, 



(2) 
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HI = 1.57 



= 2.10 



HI = 2.20 



Figure 7: The location of solutions to the non-linear equations of motion (4.21) for the 
rho meson condensate, including the new ground state, in the same units as in figure a. 
Blue curves indicate the parameter values for which A^^ {z — — oo) satisfies the boundary 
condition, while red dots do the same for = -oo). 



For the construction of the rho meson condensate, a scan is required through the space 
of C3 and coefficients in the expansion 

Again we plot the curves at which A^q\z = —00) and A^\z = —00) satisfy the boundary 
conditions, and look for intersection points. The pion condensate solution C3 = —1, = 
satisfies ( [4.21| ) for all values of this is the point on the vertical axis in figure |^. 
For m > //crit a new intersection point develops. 
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